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ABSTRACT 



An EPR study of copper halide tetrazole complexes and x- irradiated 
strontium acetate hemihydrate has been made. The powder spectra of some 
of the tetrazole complexes exhibited temperature dependent narrowing and 
broadening. The one-electron population of the copper d orbital involved 
in bonding to 2-methy 1-5-amino tetrazole has been found to be 0.56 and 
0.47, respectively for the chloride and bromide complexes. The g value 
for these two complexes are consistent with the proposed planar poly- 
meric structure for 2-substituted tetrazole-copper complexes. The 
species present in single crystals of strontium acetate hemihydrate ir- 
radiated and observed at liquid nitrogen temperature to -100° C has a 

CH 

spectrum of four lines of 1 : 3 : 3 : 1 intensity ratio and a 3 = 11.0 gauss. 
Temperature dependent exchange behavior of the spectrum for this species 
is observed and attributed to a rotating methyl group which was found to 



rotation. The species present at room temperature has an eight line 

CH 

spectrum consisting of two 1 : 3 : 3 : 1 intensity ratio patterns with a 3 = 



approximation were performed and possible fragments responsible for 



have AH^ = 3.0 + 0.3 



kca 1 
mole 




25.2 and a = 21.3 gauss. Molecular orbital calculations using the CNDO 



these spectra are suggested. 
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I . INTRODUCTION 



The work reported in this thesis consists of two parts. The first 
part, Sections II through VI inclusive, contains the results of research 
on copper halide- tetrazole complexes. The second part, which encompasses 
the rest of the thesis, examines the problem of the paramagnetic frag- 
ments produced in strontium acetate hemihydrate by x-irradiation. 

The tetrazole work was suggested by Dr. H. B. Jonassen of the 
Chemistry Department of Tulane University, New Orleans, Louisiana, to 
provide information to supplement that already gathered by other means 
to study the structure of the copper halide tetrazole complexes and the 
coordination of the tetrazoles to the copper (II) in the complex. The 
radiation damage work on strontium acetate was begun at the conclusion 
of the copper tetrazole work with the purpose of obtaining information 
about the mechanism of radiation damage in crystals. This work has 
resulted in the discovery of free radical fragments not previously re- 
ported for a crystal of this type [25]. 

For each section, the thesis is written so that the reader is given 
the theory and background necessary to understand and appreciate the 
reported results in the context of current research in the areas of 
transition metal-tetrazole complexes and radiation damage studies of 
organic crystals. 
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II. ESR OF TRANSITION METALS 



Electron spin resonance (ESR) measures the absorption spectra which 
arise from transitions between energy states. The energy states are the 
result of an interaction of the unperturbed energy states of the para- 
magnetic species with a magnetic field. In this section the theory of 
these states and their description by a spin Hamiltonian will be dis- 
cussed. The theory of the g tensor for transition metals as developed 
by Carrington and McLachlan [l] will also be considered. 

Elementary quantum mechanics provides the following relationship 
for calculating the average value of a property of a quantum mechanical 
system : 



calculated and \|f is the spatial wave function for the system. ESR is 
used to determine the energy levels of a system. An extension of the 
above relationship to a system with several energy levels can be made 
without too much difficulty. The operator for the energy is called the 
spin Hamiltonian which is, for one electron: 



In essence, the terms in the Hamiltonian describe the contribution 
to the total energy which are made by interactions of the electrons and 
nuclei with the external magnetic field and with one another. The first 
term is the electron Zeeman term. It represents the Zeeman energy, the 
energy of interaction of the electron spin moment with the external 




where 9 is an operator for the property Q whose average value is to be 



3C *0H • g • s - Sv./iH • I. + LS 

s 6 i i i 




i i 
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magnetic field. 3 is the Bohr magneton, S is the electron spin and g 
is the spectroscopic splitting tensor which may well be anisotropic. 

The second term is a summation of the Zeeman energies of the nuclei. 

1^ is the spin of the ith nucleus and is the gyromagnetic ratio of 
the ith nucleus. The third term represents the summation of the energy 
of the hyperfine interaction between the electron spin and the spin of 
the ith nucleus. is the hyperfine interaction tensor. In this 

section only the first term in the Hamiltonian will be considered: 




In the absence of spin-orbit coupling the g tensor is isotropic and 
S represents the true spin of a free electron. However, consider the 
case where g is anisotropic, since this is the case which is almost al- 
ways encountered with transition metal complexes. It naturally follows 
from above that an anisotropic g tensor implies that S cannot possibly 
represent the true spin. Instead, for anisotropic g, 

3C = 0 H • g • S 
s 

A 

where S is called the fictitious spin. 

The reader might recall the fact that in the presence of a magnetic 
field, there is an interaction between H and L, the orbital angular 
momentum of the electron. In this case, the true form of the Hamiltonian 
for the interaction of an unpaired electron with a magnetic field is 

K = 0 H • L + g p H • S 
s e 

where g is the free electron g = 2.0023. The two forms of II are 
°e s 
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equivalent. This fact is used to develop some mathematical relation- 
ships which will help in the understanding of the anisotropy of the g 
tensor . 

It can be shown that the only way an odd electron can acquire or- 

t 

bital angular momentum is through the effect of spin-orbit coupling, 
which is represented by 

-4 —> 

£ L • S 

where £ is called the spin-orbit coupling constant, the value of which 
may be extracted from atomic spectra. 

First order perturbation theory allows the calculation of modified 

—4 — ♦ 

ground state wave functions < + 1 and < - 1. £ L • S is the pertur- 

bation operator which mixes the ground state wave functions (ty o ,Q'| and 
with excited states. |cv) and |(3) are the spin functions (for 
notation, and details of development, see [l]). 

The fictitious spin can now be defined as follows: 



sj + > = %l + > sj + ) = \\ - ) 

Sj - > = -%| - ) S y | + ) = %i| - > 



In a magnetic field along the i axis (i = x,y,z) 



K = <*H. *(g. S + g. S + g. S ) . 
s 1 ix x ly y iz z 



Setting up the matrix representation of this form of K in the basis 

s 

I + ), | ~) and comparing it with the matrix representation of 

K = 0H * L + g 0H • S in the same basis, it is found that: 
s °e 
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where AE is the difference in energies between the ground state elec- 



n 



tronic level and the nth electronic level. The value of AE may be 



n 



extracted from optical spectra. 
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III. EXPERIMENTAL SPECTRA AND THEIR ANALYSIS 



A. SINGLE CRYSTAL METHOD 

Of the methods available for studying paramagnetic species in solids, 
the single crystal method provides the best results. The species to be 
studied is regularly oriented relative to the crystal axes, and this 
provides the opportunity to rotate the crystal in the magnetic field and 
measure the anisotropic interactions, giving valuable information about 
structure of the paramagnetic species. 

The type of spectrum obtained is dependent upon the orientation of 
the external magnetic field with respect to the paramagnetic species, 
and the number and type of different species in the crystal. In addi- 
tion, there may be more than one magnetically distinct site for a parti- 
cular species in the lattice. Fortunately, there is usually a parti- 
cular orientation for the magnetic field which results in a spectrum 
which is relatively easy to interpret and allows the experimentalist 
to extract the desired information. This information enables the experi- 
mentalist to find the orientation of the principal axis system of the 
species with respect to the crystal axis system. It is in the principal 
axis system that the A and g tensors are diagonal. Since each tensor 
has three diagonal components, there are three orthogonal directions for 

the magnetic field along which the values of g , g , and g may be 

xx yy z z 

extracted. A , A , and A may sometimes also be extracted, 
xx’ yy’ zz 
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B. THE POWDER AND FROZEN SOLUTION METHODS 



When single crystals of the material to be studied cannot be grown, 
or when the information cannot be extracted from single crystal studies, 
then the powder or the frozen solution "glass" method is utilized. In 
both the powder and glass, the paramagnetic species are randomly orient- 
ed with respect to the magnetic field and absorb over a wide range of 
magnetic field. The glass provides a method for diluting the species 
in order to reduce the line broadening effects which are caused by the 
interaction of neighboring complexes. 

C. LINE WIDTHS 

There are three factors which McGarvey [2] considers that can 
contribute to the width of ESR absorption lines. They are: 

(a) Anisotropy of g factors and hyperfine interactions or spin- 
interactions . 

(b) Interaction with magnetic dipoles of neighboring electronic 
and nuclear spins. 

(c) Spin-lattice relaxation times. 

The first factor is important for powder samples in which the small 
crystals are randomly oriented. In this case, the method needed to 
obtain narrower lines is to grow single crystals. If this cannot be 
done, then it is sometimes possible to extract the information from the 
shape of the broad powder spectrum. 

The second factor can often be a significant one, since magnetic 
fields produced by dipoles over distances of a few angstroms are very 
large. These large internal magnetic fields coupled with the external 
magnetic field may cause the local field at the site of the paramagnetic 
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species to vary greatly. Since the field of a magnetic dipole is in- 
versely proportional to the third power of the distance, the best way to 
reduce this broadening effect is to place the species in an isomorphic 
diamagnetic crystal or to use the glass method. 

The length of time associated with the relaxation of an excited spin 
state to a thermal equilibrium spin state given by Boltzmann statistics 
is called the spin-lattice relaxation time. The relaxation is accom- 
plished by a transfer of energy between the magnetic spins and the 
vibrational energy of the crystal lattice to which the spins are coupled. 
For systems where the spins are strongly coupled to the vibrational 
modes, the lifetime of a given magnetic state is short, resulting in 
an uncertainty in the energy which manifests itself as a broad absorp- 
tion line. This broadening mechanism is strongly dependent on the 
temperature, so that the lines broadened in this manner can be sharp- 
ened by lowering the temperature of the sample. Extremely short relax- 
ation times often occur when the paramagnetic species has an electronic 
excited state just above the ground state. In cases of this sort, it 
is often necessary to go to liquid helium temperature. 
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IV. BACKGROUND ON TRANSITION METAL- TETRAZOLE COMPLEXES 



Tetrazoles are f ive-membered , heterocyclic ring compounds which 
contain one carbon and four nitrogen atoms linked by three single and 
two double bonds. The parent compound may exist in tautomeric forms I 
and II. It was reported by Lounsburg [3] that 97% of tetrazole is of 
the form I : 




(3) 



I II 

The tetrazole ring is unusual in that the only two points of sub- 
stitution are in positions 1 or 2 and 5. Chemical investigations have 
shown that the tetrazoles are nucleophilic reagents and that the nucleo- 
philic character varies with the nature and position of the groups sub- 
stituted on the ring. 

Jarvis [4"] investigated the crystal structure of a 1:1 complex of 
copper (II) chloride and 1,2,4 triazole. He reported the unit cell 
dimension to be 

a = 6.81, b = 11.39, c = 7.13 A ; 8 = 96° 38' 

The space group is I2/c and contains four units of C 2 N 2 H 2 * CuCl^. The 
copper atom has a distorted octahedral coordination group consisting of 
two N at 1.98 A, two Cl at 2.34 A, and two Cl at 2.77 A. The structural 
unit is an infinite chain in which octahedral groups, joined by sharing 
edges, are also linked by triazole molecules. 
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H 

1,2,4 TRIAZOLE 



Jonassen [5] reported the preparation of a number of copper halide 
tetrazole complexes of differing compositions. He determined that the 
copper halide complexes of 1-substituted tetrazoles have the composi- 
tion Cu X^T^ and postulated a structure analagous to the triazole struc- 
tural unit for them. The copper halide complexes of 2-substituted 
tetrazoles have the composition CuX^T. The proposed structures of the 
1 and 2-substituted tetrazole complexes of copper (II) halides, respec- 
tively, are shown in the top and bottom of Fig. 1. Jonassen also 
carried out a reduction reaction on several of the copper tetrazole 
complexes with acetone and discovered that the product of the reaction 
was always a CuXT, regardless of the composition of the starting material. 
He reasoned that in CuXT, polymeric CuCl chains are present with bridg- 
ing tetrazole groups forming four coordinated copper (I) species. 

Garber and Brubaker [6], reporting the results of their work on 
bis (1-methy 1-5-tetrazoy 1) and bis (1-cyclohexy 1-5- tetrazoyl) nickel (II), 
stated that the reflectance spectra of these complexes indicated that 
the nickel is in an octahedral environment. They also concluded that 
the insolubility of these complexes in all common solvents suggest a 
polymeric structure. 
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Brubaker [7] carried out continuous variation experiments on 5- 
aminotetrazolato-copper (II) and showed that there is a 1:2 interaction 
of copper (II) with the tetrazole according to the mechanism: 

HT = H + + T“ K = 6.8 x 10' 7 

a 

Cu 2+ + 2T~ = CuT 2 

He determined the formation constant of the association reaction to be 
12 

K = 1.3 x 10 . Infrared analysis of this complex indicated that co- 

ordination of the tetrazole to the copper (II) does not involve the 
amino group. Jonassen [5], reporting on this same complex, stated that 
the tetrazole nitrogens form dative bonds with the copper (II) above and 
below the plane, thus forming a bridged structure. 

Bowers and Popov [8") worked on pentamethy lene tetrazole complexes 
of copper halide and compared the results with work carried out pre- 
viously on Cu(PMT),(C10 ) . They found that PMT is not a strong donor 
ligand and that at most, two PMT coordinate with Cu(II) in the halide 
complex as compared with six PMT with the perchlorate complex. In a 
few cases, they discovered solubilization of complexes when the solvent 
had strong donor properties. They found an abnormally low magnetic 
moment for the complex which is indicative of copper-copper bonding, 
but they concluded that their data did not allow an unambiguous identi- 
fication of such bonds. They also concluded, from magnetic and spectral 
evidence, that the copper (II) in Cu(PMT)X 2 is in an octahedral environ- 
ment . 

Kuska, D'ltri and Popov [9] investigated the powder spectrum of 
Cu(PMT)^ (C10^) 2 and from the observable nitrogen superhyperf ine splitting 
and comparison with optical and ESR data for other copper complexes, 
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FIGURE 1. Proposed structures of L-substituted tetrazoles (top) and 
2-substituted tetrazoles (bottom) after Jonassen [5], 
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concluded that the bonding of the PMT ligand is 80% ionic. They also 

concluded that the PMT nitrogen-copper (II) orbitals have a greater 

2 

percentage of s character than the assigned sp hybridization. 
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V. MATERIALS AND EQUIPMENT 



A. OBTAINING THE DATA 

The copper tetrazole complexes studied were suppLied by Dr. R. A. 
Henry of the Chemistry Division, Naval Weapons Center, China Lake, 
California. Their purity was checked by Galbraith Laboratories, Inc., 
Knoxville, Tennessee. The analytical data indicated that the samples 
were of sufficiently high purity for this study [5]. 

The ESR spectra were obtained with a Varian V-4502-06 X-band ESR 
spectrometer utilizing 100 kc/sec modulation with a nine inch magnet 
and field-dial regulation. The cavity was a conventional reflection 
type cavity, equipped with a Varian variable temperature control unit 
which allowed temperature regulation down to 90° K. For liquid helium 
work, the same type of spectometer modified to accomodate a liquid 
helium dewar and pumping system was used. This system allowed obser- 
vation of the sample at 3° K. 

The powder spectra were obtained by placing the powder in 4 mm 
O.D., 3 mm I.D. quartz tubing sealed at one end. The frozen solution 

% 

spectra work was accomplished in a similar manner by filling the sealed 
quartz tubing with a solution of the tetrazole complex and then freez- 
ing it with the temperature control unit. The single crystal work was 
accomplished by rotating the sample on the gear wheel of a goniometer, 
an instrument which gives the experimenter the freedom to orient the 
crystal about two axes, the axes vertical to and horizontal with the 
pole faces of the magnet. 
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The data was taken from the EPR spectra for the powder and frozen 
solution method by using the baseline of the spectrum as a reference 
and measuring the linewidths and coupling constants with a ruler gradu- 
ated in hundredths of an inch. The data for the single crystal work was 
taken from the oscilloscope display of the absorption spectrum. The 
orientation of the single crystal was changed with the goniometer while 
observing the scope display, and data taken when the absorption line 
reached an extreme or inflection point of the magnetic field at which 
resonance was occurring. The resonance fields were read off the field 
dial by centering the absorption line on the crosshairs of the oscillo- 
scope display. 

The calibration of the field dial with the linear scale of the re- 
corder was accomplished by inserting a solid mixture of 0.1% MnCl^ in 
MgO into the cavity And comparing the field dial readings and linear 
measurements of the coupling constant from the recorder with the coupl- 
ing constant as determined with an NMR probe. 

For the single crystal work, each piece of data consisted of two 
angles of rotation, a resonance field, and a linewidth. The data was 
cataloged with a stereographic net, and the three orientations at which 
data was taken for each single crystal were confirmed to be orthogonal 
to one another on the stereographic net. 

B. CATALOGING THE DATA 

Suppose a paramagnetic single crystal is to be studied. This crystal 
has a definite form which allows for definition of a three dimensional 
orthogonal crystal axis system relative to the crystal faces. The 
experimentalist has the job of finding the orientation of the principal 
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axis system with respect to the crystal axis system. Consider this pro- 

f 

cess from a different viewpoint. Reorienting the crystal in a fixed 
magnetic field is equivalent to reorienting the magnetic field with 
respect to a fixed crystal axis system. 

Consider the crystal axis system X Y Z in Fig. 2. It is easy to, see 
that all points in space may be convered by rotating a given vector 
through all angles of 0 and Suppose the experimenter wants to 

collect data from the paramagnetic crystal over several magnetic field 
orientations. The stereographic net provides a convenient method of 
cataloging the information obtained. Each value of 0 and <t> corresponds 
to one point on the stereographic net. Figure 2 gives the three dimen- 
sional picture of a principal axis system X' Y' Z' relative to the 
crystal axis system X Y Z . The representation of the same points on 
the stereographic net appears directly above for comparison. The 
stereographic net covers all points over a hemisphere. Because the 
ESR spectrum remains unchanged by a reversal of the magnetic field 
direction, corresponding points on the two hemispheres are equivalent. 

The procedure for finding X' Y' Z' consists of taking data over 
several orientations of the magnetic field and cataloging the data with 
the stereographic net. Once this is done, a reasonable guess is made 
at the three orthogonal points on the net which, in consideration of 
the data distribution on the net, seem most likely to represent X', Y', 
and Z'. The initial guesses are then fed into a computer program for 
further refinement. 

The stereographic net thus provides a convenient method for catalog- ' 
ing data and providing the experimenter with an overall picture of the 
data distribution which leads to an initial guess at the orientation of 
the principal axis system relative to the crystal axis system. 
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Wulff Net 





FIGURE 2 
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VI. EXPERIMENTAL WORK AND RESULTS 
A. FROZEN SOLUTION SPECTRA 

Frozen solution work was undertaken in order to study those com- 
plexes which did not grow single crystals or for which the powder spectra 
were unsatisfactory. The following organic solvents were tried: 

Isopentane 

Glycerine 

Benzene 

Chloroform 

Diethyl ether and mixtures of these 

Methyl ethyl ketone 

Acetone 

Ethanol 

Dimethyl formamide 

A small amount of the powder was placed in a beaker and the sol- 
vent added at room temperature. Upon very careful attempts to dis- 
solve the complexes in the solvents, it was discovered that in almost 
all attempts either the complex did not dissolve at all or what at 
first appearances seemed to be dissolution of the complex in the sol- 
vent was in actuality only a fine dispersion of the powder in the 
solvent. A few of the complexes did dissolve to a very slight degree, 
however, the spectra which were obtained in the glass at 113° K were 
impossible to interpret due to the small signal to noise ratio and 
other complicating effects. 

The only solvent which did dissolve all of the copper halide tetra- 
zole complexes was dimethylf ormamide . The copper bromide tetrazoles 
gave quite simple spectra in DMF whereas the spectra of the copper 
chloride tetrazoles showed several overlapping peaks which was indicative 
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of more than one distinct position in the frozen lattice of dimenthy Ifor- 
mamide . 

Three typical spectra of frozen DMF solutions of copper bromide 
tetrazoles are shown in Fig. 3. The top spectrum is that of copper bro- 
mide dissolved in DMF. It was used as a reference with which to compare 
the other spectra. The other two spectra exemplify, respectively, the 
phenomena of broadening of the spectral lines and shifting of the 
spectral lines. The middle spectrum is 1-benzy 1-5-amino tetrazole with 
copper (II) bromide in DMF, and the bottom spectrum is 1,5 tetramethy- 
lene tetrazole with copper (II) bromide in DMF. 

The fact that the complexes did not dissolve to an appreciable 
extent in all but DMF substantiates the conclusion of previous workers 
that the complexes have a polymeric structure which does not break down 
easily. The dissolution in DMF seems to imply that DMF coordinates 
more strongly to copper than do the tetrazoles. The data in Table 1 
leads to a partial confirmation of this fact, and also allows for some 
other conclusions. 

Previous workers [9^ concluded that the bonding of the PMT ligand 
to copper is 80% ionic and that PMT is a poor donor. Data in Table 1 
indicates that the PMT ligand has been replaced by DMF and that the 
same thing has happened for some of the other complexes. This is not, 
however, the case for all the complexes as can be seen. For the two 
complexes which exhibit values of AH and A different from those listed 
for the copper bromide case, it must be concluded that the tetrazole 
ligands have not been completely replaced by the DMF, but rather, a 
type of coordination still exists between the copper and the tetrazoles. 
The difference in AH values is due to the difference in the ligand 
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COMPLEX 



CuBr^ 

2 phenyl 

1 Methy 1-5-ethyl 
1,5 pentamethylene 

2 methy 1-5-amino 
1 benzy 1-5-amino 
1 methy 1-5-amino 
1,5 tetramethylene 
1,5 trimethylene 



AH (gauss) 
448 + 6 
442 + 6 
448 + 6 
448 + 8 



AVERAGE A (gauss) * 
11.9 + 6 3 
119 + 6 3 
116 + 6 3 
119 + 8 2 



lines too broad to be observed 



lines too broad to be observed 
lines too broad to be observed 
416 + 8 
412 + 8 



139 + 8 2 
138 + 8 2 



’’^Superscript indicates the number of A values which were averaged. 

In cases where less than three values were averaged, the remaining 
A values were lost to line broadening effects. 




TABLE 1. CuBr -Tetrazole Complexes in Frozen Solution of Dimethylfor- 
mami3e 
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splitting of the electronic levels of the copper ion. Since the AH 
values are inversely proportional to the ligand splitting, AE (see 
Introduction), larger AH values imply a smaller splitting by the ligand. 
The A values are a measure of the localization of the electron in the 
copper orbital. The larger the A value, the more the electron is 
localized. 

B . POWDER SPECTRA 

Due to the lack of success in the attempts to grow single crystals 
of the complexes supplied, spectra were taken of the complexes in 
their powder or microcrystalline form. The resulting spectra all 
exhibited extremely broad lines at room temperature. Initially attri- 
buting the broad lines simply to spin- lattice relaxation effects, a 
series of experiments was performed on four particular complexes which 
exhibited anisotropic line shapes at room temperature. These experi- 
ments revealed temperature dependent behavior of the spectra of all 
four complexes. Figure 4 shows the variation of the powder spectra of 
the four different tetrazole complexes with temperature. The four com- 
plexes studied were: 

Figure 4A: 1-methy 1-5-ethyl tetrazole with CuCl^ 

Figure 4B: 1-methy 1-5-amino tetrazole with CuCl^ 

Figure 4C: 1,5 -pen tame thy lene tetrazole with CuBr^ 

Figure 4D: 1-benzy 1-5-amino tetrazole with Cu Br^ 

Encouraged by the fact that for some of the complexes, the spectral 
lines narrowed at lower temperature, it was decided to try work at 
liquid helium temperature. Initial preparations with a spectrometer 
similar to that used for the liquid nitrogen work proceeded slowly and 
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373 ° 
273 ° 




figure: 4 a 
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FIGURE"- 4 B 
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FIGURE 4-0 
373 ° 
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with much difficulty, but finally a powder spectrum of 1-methy 1-5-amino 
tetrazole copper (II) chloride complex was obtained at liquid helium 
temperature. The result is shown in Fig. 4A. It was concluded on the 
basis of this experiment that the liquid helium work would bear little 
fruitful results, and so efforts in this direction were terminated. 

Room temperature powder spectra of the complexes which exhibited 
isotropic line shapes were taken and the line width measured. Room 
temperature powder spectra of the three copper (I) complexes were also 
taken. The results are listed in Table 2. 

Although no positive explanation can be given for the temperature 
dependent behavior of some of the complexes, hypotheses may be advanced 
to explain the effects. Typical behavior is exemplified by l-methyl-5- 
ethyl tetrazole with CuCl^ in Fig. 4A. The spectrum seems to alter- 
nately broaden and narrow as the temperature is raised or lowered. The 
importance of some line broadening effects will now be considered to 
explain this behavior. 

At high temperature, the effects of spin lattice relaxation are 
important. In addition, there exists the possibility of rapid exchange 
of the tetrazole ligand between two positions it could occupy (Fig. 1, 
top) in the puckered polymeric chain. This possibility is advanced in 
light of the discovery of previous workers [9] that the tetrazole ligand 
is in general a poor donor, and therefore does not interact strongly 
with the central copper. The polymeric chain would then be free to 
dynamically "repucker" by making and breaking bonds with the tetrazole. 
Spin- lattice effects would tend to broaden the spectral line and ex- 
change tend to narrow it. At any particular temperature the dominance 
of one effect or the other would depend upon the structure of the 
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TABLE 2. Powder Spectra of Copper Halide Tetrazole Complexes 



complex and the temperature dependence of the exchange mechanism. At 
lower temperatures, where the spin- lattice effects are greatly diminished 
and the lattice vibrations are very small, it is possible for the copper 
atoms to be close enough so that the electron is affected by both copper 
nuclei. This effect leads to dipolar relaxation and line broadening. 

C. SINGLE CRYSTAL WORK 

Of the many copper halide tetrazole complexes supplied by R. A. 

Henry of China Lake, only the 2-methy 1-5-amino copper (II) chloride and 
bromide tetrazole complexes were in the form of crystalline chips large 
enough to be treated as single crystals. Numerous attempts at growing 
single crystals of the other canplexes were made, and all but one at- 
tempt met with failure. The conventional method of growing crystals by 
slow evaporation from solution was impractical because the copper 
halide tetrazole complexes which were supplied were insoluble in all 
common organic solvents tried. Basically, all of the techniques for 
growing single crystals were based on the method of slow diffusion of 
the copper halide and the tetrazole dissolved in a suitable solvent. 

The slow diffusion technique was necessary because of the tendency of 
the complexes to immediately precipitate when a solution of tetrazole in 
solvent was poured into a solution of copper halide in the same solvent. 
The complexes have a very large formation constant [7], 

The successful attempt at growing a single crystal of 2-methyl-5- 
amino copper chloride tetrazole complex yielded a crystal approximately 
ten times the size of any of the small crystals initially supplied from 
China Lake. A few grams of copper chloride was placed in the bottom of 
a long pyrex tube sealed at one end. The tube was filled with glass 
wool and a few grams of the 2-methy 1-5-amino tetrazole placed on top of 
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the column of glass wool. The entire tube was filled with absolute 
ethanol and then stoppered and allowed to stand at room temperature. 
After a period of two months, a crystal about the size of a pea was ob- 
served to be clinging to the column of glass wool. Due to the time 
period involved in this technique, it was not tried with other tetra- 
zoles. Of all the techniques for growing the single crystals of copper 
halide tetrazoles, this one showed the most promise. However, the com- 
plex from which the crystal was grown had already grown smaller single 
crystals. This was not the case for the rest of the complexes which 
came in powder or microcrystalline form. 

The results of the single crystal work with the two complexes men- 
tioned above is reported in Table 3. The spectra of the single crystals 
showed one sharp, isotropic line at any orientation of the magnetic 
field which means that there exists only one magnetically distinct site 
in the crystal. In both crystals, the g values are all different with 
two out of each set differing only by 0.03. Both of these observations 
are consistent with the proposed polymeric structure of a 2-subs tituted 
tetrazole (Fig. 1, bottom). The two g values which are almost equal 
correspond to the g values taken when the magnetic field is parallel 
to or perpendicular with the plane of the polymer. The proposed struct- 
ure also has only one magnetically distinct site. 

The g values may be used to make a rough calculation as to the 
amount of covalency in the copper- tetrazole bond. Assuming the struct- 
ure at the bottom of Fig. 1 to be the structure of the complexes under 
consideration here, two formulae for the g values of an octahedral 
copper complex with a tetragonal distortion corresponding to elonga- 
tion along a four-fold symmetry axis may be developed. The development 
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is similar to the development of the g value formuLae derived in the 
introduction. The final formulae are: 



g = 2.002 + P 

g = 2.002 + P 



5 = the spin orbit constant for 
copper = 829 cm"^ 

P = the one electron popula- 
tion of the copper d 
orbital [ 2 ] 



Jonassen [5] reported the values for for 
value for A^. With the above formulae, P may 
used to predict A^. The results are given in 



both complexes but 
be calculated and 
Table 3. 



gave no 
then 
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Position and kind 
of substituent on 
tetrazole 


Composition 


g values 


linewidths 

(gauss) 


2-methyl-5 amino 


2(C 2 H 5 N 5 ) CuC1 2 


2.043 + 0.002 


24.5 + 0.5 






2.073 


29.4 






2.266 


49.0 


2-methy 1-5 amino 


2(C 2 H 5 N 5 } CuBr 2 


2.043 + 0.002 


22.0 + 0.5 






2.073 


34.3 






2.218 


36.8 



Now use 




- 2 . 002) A^ 

Ik 



( 1 ) 



In the computation of 



a 2SP 

2 (g L - 2.002) 



( 2 ) 



Chloride Complex 
Given: = 2.266 

A = 14100 cm" 1 

Obtain from (l) P = 0.56 



Bromide Complex 
8 U = 2.218 

A = 14380 cm" 

P = 0.47 



Using these values of P and taking g = 



average of g and g 

xx yy 



= 2.058 



Result: A^ = 16580 cm 



A 2 = 13920 cm 



TABLE 3. RESULTS OF SINGLE CRYSTAL WORK 
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VII. THEORY OF THE ESR OF TRAPPED ORGANIC RADICALS 



The effective spin Hamiltonian describing the interaction of the 
nuclei and electron of a trapped organic radical with the magnetic 
field and with one another has previously been written down in the 
section on the ESR of transition metals. It will be listed again to 
aid in the following discussion: 



3C = 3 H • g • S - Ev.SH - I. + LS - A • I. 

i 1 1 , i 1 

In the section on transition metals, the first term of the effective 
spin Hamiltonian was thoroughly described. For trapped organic radi- 
cals, the g tensor of this term deviates only slightly from the free 
spin g value. It is symmetrical and its principal values are always 
close to the free spin g value. The second term represents the inter- 
action of the nuclear spins with the magnetic field. The effect of 

this term is primarily to relax the selection rule Am^. = 0. 

* 

Reference 10 contains an excellent description of the third term 
which represents the hyperfine interaction between the electron spin S 
and the nuclear spin I described by the second-rank tensor A. Insofar 
as g is nearly isotropic, the hyperfine tensor A will be a symmetrical 
3X3 matrix which can be diagonalized and resolved into isotropic and 
anisotropic components. The effective Hamiltonian H' for the hyper- 
fine interaction energy can be written: 



t|r 2 (0> S • I - gpy 



- . - -3 

S Ir 



3C* = 

The eigenvalues of H' are the principal values of A 
and anisotropic components of A correlate with the 
terms, respectively. 



- 3(S • r)(I • r)r 5 
, and the isotropic 
first and second 
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The first term in the equation for H' gives rise to a nonzero hyper- 
fine splitting only if there is a finite probability of finding the un- 
paired electron at the interacting nucleus, as indicated by the factor 
2 

ijf (0). For this reason the isotropic interaction is usually called the 
contact or Fermi interaction. It is apparent that the isotropic hyper- 
fine interaction A will only be nonzero if \lr (o) is nonzero; that is, if 
the orbital of the unpaired electron has s-character. 

The second term in the equation for H' imparts to the effective 
Hamiltonian an anisotropic or direction dependent character. This 
anisotropic part of the hyperfine tensor is a traceless tensor which 
describes the hyperfine interaction energy consequent upon the classi- 
cal interaction of the electronic and nuclear dipoles. 

In general, the hyperfine interaction of n nuclei of spin 1/2 will 
give rise to a spectrum consisting of 2° equally intense lines. If two 
or more equivalent spin 1/2 nuclei interact, superposition of some of 
the lines results in a spectrum of n + 1 lines of binomial intensity 
distribution. 
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VIII. BACKGROUND ON THE ESR SPECTRA OF TRAPPED ORGANIC RADICALS 



A considerable amount of work has been done on the subject of the 
ESR spectra of radiation damaged organic crystals. Naturally, the work 
has covered those materials from which single crystals are easily ob- 
tained. Crystals of saturated carboxylic acids and their salts, the 
simpler amino acids and certain unsaturated acids, when irradiated with 
x-rays, gamma rays, or high energy electrons, have been found to pro- 
vide suitable hosts for the trapped radicals produced by such irradia- 
tion. Analysis of the spectra from radiation damaged crystals of these 
materials is primarily concerned with the identification of the trapped 
radical(s) via the interpretation of the proton hyperfine interaction. 
One type of radical which is commonly detected is a n - electron radi- 
cal centered on a carbon atom. In such a radical the unpaired elect- 
ron predominantly occupies a carbon 2p orbital perpendicular to the 
i - 

plane of the radical. 

Kispert and Rogers [13] reported the interpretation of the ESR 
spectrum of single crystals of sodium acetate trihydrate irradiated 
with 1 Mev electrons at 77° K. At this temperature, they observed a 
four line pattern with intensity ratio of 1 : 3 : 3 : 1 and proton hyper- 
fine splitting of 22.5 gauss. Each of these lines had two satellites 
spaced 5.6 gauss on either side to form a twelve line pattern. They 
interpreted the main four line pattern to be due to the methyl radical, 
•CH^. The satellite lines were attributed to spin flip transitions 
arising from a simultaneous change in the spin state of the odd elect- 
ron on the methyl radical and of the neighboring protons of the 
molecules of water of crystallization. One of the first workers to 
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investigate the ESR of an x-irradiated single crystal of methyl malonic 
acid was Heller [14], He irradiated the malonic acid at room tempera- 
ture and observed the spectra at room temperature, 77°K, and 4° K. 

Heller interpreted his spectra in terms of the two species CH^CCCOOH)^ 
and CH^CHCOOH. He concentrated his efforts on the first and con- 
cluded, from the spectra, that the three methyl protons are equivalent, 
and the principal values of their hyperfine coupling tensor are 75.4 
(along the C-CH^ bond), 68.8 and 68.6. He also concluded that the 
isotropic component of the coupling tensor (and, therefore, the spin 
density on the methyl protons) is positive, and that the methyl group 
executes nearly free rotation about the C-C bond at 4.2° K. Tamura, 
Collins, and Whiffen [15] irradiated various malonic acids at 77° K 
and observed the ESR spectra as the temperature was raised. They 
found that the radical formed at 77° K was of the form R-CH2 _ CH(C00H) , 

and that as the temperature was raised, a reaction took place in which 

# 

the original radical disappeared and a radical of the form RC^C^OOH)^ 
showed up in its place. 

J. R. Morton [l6"j reported the ESR spectrum of (CH^^CO^ trapped 
in O'-aminoisobutyric acid at 40° K. He concluded that at this tempera- 
ture, one of the methyl groups is undergoing virtually free rotation 
about the C-CH^ bond, while the other methyl group is undergoing hinder- 
ed rotation. Box and Freund [17] investigated x-irradiated cv-amino- 
isobutyric acid irradiated at 77°K. They allowed the crystal to warm 
gradually and observed the resulting spectra. They detected distinct 
intermediate ESR spectra due to preliminary conformations of the radi- 
cal present at room temperature. They found the intermediate spectra 
to be better resolved and easier to interpret than the final spectrum. 
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They also found that cooling the intermediate fragment to 77° K showed 
no change in the spectrum, but that cooling the room temperature frag- 
ment did show a changing spectrum. They interpreted this behavior as 
resulting from different conformations of (CH^^CCOOH. 

Miyagawa and Gordy [18] first reported the results of their investi 

gation of. room temperature irradiated single crystals of alanine. They 

proved that the stable free radical produced by the irradiation was of 

the form CH^CHR, where R is a group which has no nuclei with detectable 

coupling. The hydrogens of the methyl group of the radical were shown 

to have equivalent, isotropic coupling of 26 gauss each, essentially 

independent of the frequency of observation. Horsfield, Morton, and 

Whiffen [19 and 20] observed the change in the ESR spectrum of the room 

temperature gamma irradiated single crystal of l-a-alanine on cooling 
o 

to 77 K, and interpreted their results in terms of the cessation of 
rotation of the methyl group of the trapped radical CH^CHCOOH. The 
isotropic hyperfine couplings of the three hydrogens of the methyl 
group were found to be 120, 76, and 14 Mc/sec compared to the value of 
70 Mc/sec each at 300° K. In the intermediate temperature range bet- 
ween 100° K and 200° K, the spectra were found to differ from the 
spectra above and below these temperatures. From the line width change 
near 100° K and 200° K, an estimate of approximately 4 kcals/mole was 
made for the activation energy hindering the rotation of the methyl 
group. Miyagawa and Itoh [21] determined that room temperature gamma 
irradiated single crystlas of dl-alanine grown from heavy water ex- 
hibited the same hyperfine structure as those crystals grown from nor- 
mal water except that the line widths were considerably sharper. They 
observed the spectra as the temperature was lowered and concluded that 
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the rotation of the methyl group was hindered at lower temperatures. 

They used a modified Bloch equation to estimate the frequency of rota- 
tion of the methyl group from the absorption patterns at various temper- 
atures, and by applying the classical theory of rate processes to 
analyse the relation between the frequency of rotation and the tempera- 
ture, they found the potential barrier to internal rotation of the methyl 
group to be 3.6 kcal/mole. 

Sinclair and Hanna [22 and 23] examined the ESR spectra of 1-alanine- 
13 o 

1- C, x-irradiated at 80 K. They concluded that the low temperature 

initial species has a high unpaired electron population on the carboxyl 

carbon and that as the low temperature radical is allowed to decay 

thermally to the fragment CH^CHCOOH at 140-150° K, there is a transfer 

of the unpaired electron from the carboxyl group to the Q’-carbon. The 

CH^CHCOOH is at a different orientation in the crystal than the 

CH^CHCOOH formed by irradiation at room temperature. Upon warming to 

room temperature, a reorientation of the fragment occured which gave 

the usual room temperature ESR spectrum. A feature of the intermediate 

fragment was that the methyl protons remained equivalent even when re- 

cooled to 80° K. This was in contrast to the room temperature fragment 

where the rapid interconversion of methyl protons was frozen out by 
o 

cooling to 77 K. 

Wells and Box [24] studied the rotation of the methyl group in the 
free radical CjeCH 2 -CONH-C(CH 3 >COOH by ESR and ENDOR spectroscopy. They 
analyzed their results using a quantum mechanical hindered rotor model 
with Mathieu functions as basis functions. Their experimental results 
indicated a barrier height of less than 0.5 kcal/mole for the methyl 
group rotation. 
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In recent paper on x-irradiated single crystals of zinc acetate 
dihydrate, Tolies, Crawford, and Valenti [25] reported the presence 
of the methyl radical in the case of irradiation at 77 1 K, with the 
subsequent disappearance of this radical and the appearance of 
as the temperature was raised. They proposed a mechanism to explain 
the decay of the first species and corresponding increase in the 
second species. The mechanism involves abstraction of a hydrogen atom 
by the methyl radical from a neighboring acetate ion. They supported 



this hypothesis by synthesizing Zn(CH o DC0 o ) • 2H o 0 , , , ^ . 

• r ° 2 2 2 2 and showed that the 

decay and growth of the above species exhibits an isotope effect. 

Spectra taken over a range of temperatures revealed both the slow and 

rapid exchange of the hydrogen atoms seen previously [e.g. 16,17,18], 

with the barrier to internal rotation of the group being 5 kcal/ 

mole. / 
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IX . EXPERIMENTAL WORK 



A. GROWING OF CRYSTALS 

Single crystals of strontium acetate hemihydrate were grown by slow 
evaporation from a saturated aqueous solution at room temperature. 
Crystals in which heavy water was exchanged for the water of hydration 
were grown by making a saturated solution of the strontium acetate 
hemihydrate in 99 % D^O at room temperature, and then allowing slow 
evaporation of the solution. The single crystals crystallized in the 
form shown in Fig. 5 which enabled definition of an orthogonal crystal 
axis system. The crystals cleave easily parallel to the XY plane. 

The XYZ crystal axis system is shown in Fig. 5. The numbers along the 
edges of the crystal are the angles the crystal faces make with the 
XY plane. 

B. IRRADIATION OF THE CRYSTAL 

The single crystals were exposed to x- irradiat ion at 50 KV and 40 
ma from a copper target. The crystals irradiated at room temperature 
were placed on scotch tape and mounted in the open windows of the x-ray 
machine allowing the port covers to be closed at all times. For room 
temperature irradiations, the crystals were exposed for one hour. For 
irradiation at liquid nitrogen temperature, the crystals were placed in 
a quartz tube sealed at one end. The other end was stoppered to pre- 
vent condensation of liquid air into the tube. The tube with crystals 
was placed in a styrofoam ice bucket and precooled gradually by slowly 
pouring liquid nitrogen into the bucket and allowing the cold gas to 
flow over the tube. Precooling was necessary to prevent shattering of 
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FIGURE 5. 



A single crystal of Sr(CH o C0.) • %H 0 showing the angles 

J A 2 ^ 

between the edges and faces of the crystal, and the ortho- 
gonal crystal axis system. 
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the crystal when exposed to a sudden temperature change. Eventually, 
the tube was immersed completely in liquid nitrogen. The bucket was 
mounted on a ring stand with a lead shield attached, and a wooden cover 
with a hole drilled in the top of the bucket. The bucket was placed 
in front of the x-ray port, and the tube positioned in the x-ray beam 
so that the crystal received maximum dosage for a two hour period. The 
position of the tube in the beam was checked by darkening the room and 
lowering a piece of fluorescent paper into the bucket. The tube and 
crystal formed a shadow on the fluorescent paper. 

C. CRYSTAL ALIGNMENT 

The orientation of the crystal in the spectrometer cavity had to 
be determined as precisely as possible in order to determine the loca- 
tion of the magnetic field with respect to the crystal axis system. 

The crystal was mounted on a quartz rod which had two flats ground 
parallel to each other on both ends of the rod. The crystal was held 
on the rod by means of a platinum spring attached to the rod by a 
piece of scotch tape. The end opposite the mounted crystal was clip- 
ped to an angle indicator positioned over the center of the spectro- 
meter cavity. The crystal was positioned on the rod by use of a micro- 
scope, one of the eyepieces of which had crosshairs. The orientations 
of the crystal on the rod were chosen so that the well defined edges of 
crystal were parallel or perpendicular to the edges of the rod. The 
angles between the faces and edges of the crystal were determined with 
an optical goniometer. For liquid nitrogen work, the crystal alignment 
was accomplished without warmup of the crystal by keeping the crystal 
immersed in a bath of liquid nitrogen under the microscope. 
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D. TAKING THE SPECTRA 



The spectra were taken with the spectrometer previously described 
for the Liquid nitrogen work on copper tetrazole. The temperature was 
controlled with the temperature control unit attachment to the spectro- 
meter used for the copper tetrazole work. 
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X. RESULTS 



A. ROOM TEMPERATURE SPECIES 

Spectra obtained from x-ray damaged single crystals of strontium 
acetate hemihydrate at room temperature were found to depend upon the 
history of the crystal as well as the orientation of the magnetic field 
relative to the crystal axis system. As can be seen in Fig. 6 (top), 
crystals irradiated at liquid nitrogen temperature, slowly warmed to 
room temperature, and then observed at room temperature did not exhibit 
the extra central line in the spectrum as did the crystals irradiated 
and observed at room temperature (Fig. 6, bottom). Both spectra were 
obtained at the same crystal orientation of <f> = 0°, 0 = 90°. Reorient- 
ation of the magnetic field resulted in a coalescence of the two 
1:3:3:1 patterns into what appears to be a 1:4:6:4:1 pattern as in 
Fig. 7. The orientations for the magnetic field are from top to bottom, 
respectively ^ = 0°, 0 = 45°; <t> = 0°, 0 = 80°; <t> = 0°, 0 = 90°. The 
spectra can be interpreted in terms of a free radical species with 
three equivalent hydrogens split by a fourth, inequivalent hydrogen. 

The spectra obtained also indicate that there is only one magnetically 
distinct site in the crystal for the room temperature species. 

For some orientations of the magnetic field, the base line of the 
spectrum was not straight, but rather convex upward. The sweep range 
was therefore changed from 250 gauss to 1000 gauss and the spectra in 
Fig. 8 were obtained. The spectra indicate the possible presence of 
a triplet. 
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FIGURE 6 
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FIGURE 8. Broad line species seen in SrCCH^COp^ * ^ H 2 ° irradiated and 
observed at room temperature. 
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Several spectra of the room temperature species were taken at 
different magnetic field orientations. The coupling constant for the 
inequivalent hydrogen was not always obtainable due to the overlapping 
of the spectral lines at several orientations. However, the coupling 
constant for the CH^ group was always obtainable. The methyl group 
coupling constants and corresponding magnetic field orientations were 
plotted on a stereographic net, and an estimate of the positions of 
the principal axis system and the coupling constants along the princi- 
pal axis was made, assuming axial symmetry for the rotating methyl 
group. The results are: 

= 27.1 + 0.4 gauss, direction cosines (-0.32, 0.38,-0.87) 

A = 24.2 + 0.4 gauss 
J. 

The errors in the direction of the principal axis are estimated to be 
ten degrees and for the coupling constants, 0.5 gauss. The coupling 
constants and magnetic field orientations for the inequivalent hydrogen 
were least squared [26]. The results are, for the three principal 
values of the A tensor: 

(a) 11.61 + 0.40 gauss 

(b) 32.98 

(c) 19.32 

The direction cosines for the principal axes are: 

(a) (0.72, 0.69, 0.80) 

(b) (-0.18, 0.73, 0.98) 

(c) (0.67, -0.72, 0.18) 

in the XYZ axis system in Fig. 5. The data used is in Table 4. Figure 
9 shows the angular dependence of the coupling constant. The graphs 
indicate that the methyl group coupling constants show a cos9 dependence. 
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0 25 

30 

35 

40 

45 

50 

55 

60 

65 

70 

80 



a 3 (gauss) 

25.54 + 0.35 

25.24 

25.24 

25.06 

24.90 

24.56 

24.56 

24.38 

24.38 

24.20 

24.38 



70 30 25.58 

45 24 . 90 

60 24.38 

75 24.38 

90 24.20 

105 24.90 

120 25.41 

135 26.28 

150 26.45 



45 0 26.10 

15 25.93 
30 25.54 
45 25.24 



TABLE 4. Coupling constants and magnetic field orientations for the 
room temperature species. 
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Table 4 continued 



j 

(45) 



90 



25 



0 

60 

75 

90 

105 

120 

135 

150 



a 3 (gauss) 
24.55 + 0.35 
24.38 
24.20 
24.20 
25.24 
25.76 
26.28 



0 

15 

30 

45 

60 

75 

90 

105 

120 

135 

150 

165 



26.28 

25.93 

25.24 

25.06 

24.90 

24.38 

24.55 

25.06 

25.59 

26.10 

26.45 

26.23 



30 

45 

60 

75 

90 

105 

120 

135 

150 

165 



25.42 

24.72 

24.38 

24.38 

24.20 

25.06 

25.42 

26.28 

26.45 

26.80 
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Table 4 continued. 



j 

115 



9 a 3 (gauss) 



0 


25.93+0.35 


15 


25.93 


30 


25.54 


45 


25.06 


60 


24.38 


75 


24.77 


90 


24.56 


105 


24.90 


120 


25.93 


135 


26.45 


150 


26.97 


165 


26.62 



e 


a (gauss) 


25 


17.12 + 0.35 


30 


15.04 


35 


13.83 


40 


12.10 


45 


11.24 


50 


11.06 


55 


11.93 


60 


13.83 


65 


15.73 


70 


16.94 


30 


L5.22 


45 


L0.54 


60 


12.44 


75 


L7.29 


120 


29.39 


135 


30.60 


150 


30.26 
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Table 4 continued. 





e 


(gauss) 


45 


15 


20.74 + 0.35 




30 


16.94 




45 


13.31 




60 


12.80 




75 


15.90 




90 


20.22 


70 


30 


20.22 




45 


18.15 




60 


16.94 




75 


17.46 




90 


14.01 
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a in gauss 
H 6 




FIGURE 9A. H atom coupling constant versus orientation of magnetic 
field in XYZ axis system. 0=0° 
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a in gauss 
H 




FIGURE 9B . <i> = 25° 
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a in gauss 
H 




FIGURE 9C . p = 45° 



62 



a^ in gauss 




FIGURE 9D. = 70° 
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in gauss 



CH, 




FIGURE 9E . 



CH^ group coupling constant versus orientation of 
magnetic field in XYZ axis system. </> = 0° 
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a 



CH 



3 



in gauss 




FIGURE 9F. </> = 25° 



65 



in gauss 



CH, 




FIGURE 9G. 4 - 45 
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in gauss 





FIGURE 9H . <f> = 70° 
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in gauss 



a 



CH 



3 




FIGURE 91. ^ = 90° 
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in gauss 



a 








— 23.5 



i 



120 90 60 

J I L 

9 in degrees 



1 



i 



FIGURE 9J. i = 115° 
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Powdered strontium acetate hemihydrate was recrystallized from 997> 
D^O and the resulting radiation damaged crystals showed the same spect- 
rum as single crystals recrystallized from normal water. It may be con 
eluded that the water of hydration does not enter into reactions involv 
ing the initial free radical fragments. 

An interesting sidelight to the room temperature work is that a 
crystal of zinc acetate dihydrate which had been x-ray damaged and then 
set aside for one month was looked at on the spectrometer. Surpris- 
ingly, the spectrum obtained was the same 1:4: 6:4:1 pattern obtained 
for strontium acetate hemihydrate. Apparently, the radical species 
identified by Tolies, Crawford, and Valenti [25] was an intermediate 
in a reaction with a very small rate constant or large activation 
energy which leads to a final product similar to that in strontium 
acetate. The spectrum of the species in zinc acetate dihydrate is 
shown in Fig. 10. 

B. INTERMEDIATE SPECIES 

Single crystals of strontium acetate hemihydrate irradiated at 
liquid nitrogen temperature and observed at -100° C exhibited a 1 : 3 : 3 : 1 
pattern which split into eight lines of equal intensity as the tempera- 
ture was lowered to -160° C. Spectra for the species present at -100°C 
were taken over several orientations of the magnetic field, the coupling 
constants and magnetic field orientations plotted on a stereographic 
net, and, assuming axial symmetry for a rotating methyl group, the 
hyperfine coupling constants and direction cosines of the species rela- 
tive to the XYZ axis system in Fig. 5 were calculated to be [26]: 
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A = 12.84 + 0.30 gauss 

A^ = 10.04 + 0.30 gauss 

X = (0.81, 0.58, 0.06) 

Y = (-0.59, 0.80, 0.12) 

Z = (0.02, -0.14, 0.99) 

The data fed into the program is in Table 5. The average isotropic 
hyperfine coupling constant is 9.2 gauss. 



71 



FIGURE 10. Spectrum of one month old x-ray damaged single crystal of Zn(CH 3 C0 2 ) 2 * 2H 0 
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i. 9 a 3(Mc/sec) 



0 


35.89 + 1.00 


20 


36.83 


40 


36.27 


60 


32.16 


70 


31.78 


80 


30.28 


90 


29.35 


100 


27.86 


110 


28.60 


120 


28.42 


140 


31.03 


160 


33.84 


0 


35.33 


20 


36.36 


40 


34.96 


60 


31.03 


70 


29.16 


80 


28.42 


90 


27.86 


100 


27.29 


110 


28.04 


120 


28.98 


140 


31.22 


160 


34.21 



TABLE 5. Coupling constants and magnetic field orientations for the 
-100° C species. 
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Table 5 continued 



l e_ 

45 0 

20 

40 

60 

80 

90 

100 

110 

120 

140 

160 



a 3 (Mc/sec) 

35.71 + 0.35 

34.40 

31.41 
28.79 
27.48 
27.11 
28.60 

29.72 
31.59 
33.46 
35.33 



135 0 35.71 

20 36.64 

40 35.33 

60 33.28 

70 30.10 

80 29.54 

90 28.23 

100 28.04 

110 27.86 

120 28.60 

140 28.60 

160 32.90 
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The spectra of this same species were obtained as a function of 
temperature. The computer program following the appendix, which is 



based upon the mathematical model for motional and exchange narrowing 
of spectral lines, was used to predict the spectrum of the rotating 
methyl group for various rotational frequencies. In the drawing be- 
low, the eight line pattern resulting from splitting by three equiva- 
lent protons corresponds to the zero rotation case. In the case of 
rotation of the methyl group, all of the lines except the two outside 
lines exchange, and the extent to which they exchange is a function of 
the frequency at which the methyl group is rotating: 



Various rotational frequencies were fed into the above mentioned com- 
puter program, and the computed spectra compared with the experimental 
spectra. It was possible to obtain the frequency of rotation as a 
function of temperature from this method, and an Arrhenius plot enables 
the determination of the AH^ and AS^ for internal rotation of the methyl 
group. The results are: 




13313 



aaa 



aw(3 a$a 



aBB Ba'Q' 



Bek'S BBCk 1 



AH^ = 3.0 + 0.3 kcal/mole 
AS 5 ^ = 6.0 + 0.3 cal/mole deg 
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The predicted and experimental spectra for a constant magnetic field 
orientation and varying temperature are shown in Fig. 11A. The dotted 
line spectra are the calculated spectra. Figure 11B is the Arrhenius 
plot. . ; 

C. LOW TEMPERATURE OBSERVATION 

Spectra taken for single crystals of strontium acetate hemihydrate 

irradiated and observed at liquid nitrogen temperature without warmup 

show a multiline pattern, as can be seen in Fig. 12, which does not 

o 

show up again after the crystal has been warmed to -100 C and then re- 
cooled to -160° C. The orientation of the magnetic field is <t> - 45°, 

9 = 70°. 

D. MOLECULAR ORBITAL CALCULATIONS 

Molecular orbital calculations for CH^CO^, CH^CO, and CH^CHO based 
on the CNDO approximation were made [27,28]. For CH^CO^, the calcula- 
tion predicts the unpaired spin density to be localized in the two 
oxygen 2p orbitals, with 0.499 spin density in each. For CH-CO, the 

Z j 

calculation predicts an average spin density on the methyl carbons of 
0.01916 which, when multiplied by 500 gauss which is the hydrogen atom 
coupling constant, yields a 9.6 gauss coupling constant. The spin den- 
sities are, for the carboxyl carbon, 0.57; on the oxygen, 0.32 and on 
the methyl carbon, 0.09. The calculation for CH^CHO reveals that the 
average spin density on the methyl hydrogens is 0.014076 which gives a 
coupling constant of 23.4 gauss. The spin density on the oxygen is 0.35 
and on the carboxyl carbon, 0.47. The isotropic coupling constant 
calculated from McConnell's relation [l] for the inequivalent hydrogen 
is 10.7 gauss. 
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FIGURE 11 B. Arrhenius Plot for determing activation Energy 
for rotating methyl group in -100° C species. 
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FIGURE 12 
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E. SUMMARY AND CONCLUSIONS 



Radiation damage of strontium acetate hemihydrate single crystals 

initiates a multistep solid state free radical chemical reaction. At 

liquid nitrogen temperature to -100° C, a radical of the form CH^ - X 
CH 

with a 3 = 9.2 gauss and a pattern of eight equally intense, equally 

spaced lines exists. The lines exhibit exchange effects to become at 

-100° C a four line 1 : 3 : 3 : 1 pattern. The exchange is attributed to a 

rotating methyl group with barrier to internal rotation equal to 3.0 

kcal/mole. At room temperature, the species is of the form CH^CH-X 
CH H 

with a 3 = about 25.2 gauss and a = about 20. gauss. The spectrum 
of the room temperature species is two 1 : 3 : 3 : 1 patterns offset from 
one another because of splitting by a hydrogen atom. In addition, a 
broad line species exists at room temperature which may be a triplet 
or a strontium ion. 

Comparison of the above results with the molecular qrbital calcula- 
tions shows that the methyl group coupling constant of 9.6 gauss calcu- 
lated for CH^CO compares favorably with the 9.2 gauss methyl group 
coupling constant determined experimentally for the low temperature 
species. CH^CO would also be expected to give the same spectrum as 
the low temperature species. The calculated methyl group coupling con- 
stant of 23.4 gauss for CH^CHO compares favorably with the approximate 
25.2 gauss coupling constant for the room temperature species. CH^CHO 
would give the pattern described for the room temperature species. 

CH^CO^ is ruled out as a possibility because of its spin density distri- 

* 

bution on the two oxygens. Confirmation of the identity of the species 

in the crystal might be made by irradiating a single crystal of 
13 

SrCCH^ CC^) 2*^2© anc * su bsequent ly examining the ESR spectrum. 
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APPENDIX A 



A MATHEMATICAL MODEL FOR THE MOTIONAL 
AND EXCHANGE NARROWING OF SPECTRAL LINES 

The temperature dependent behavior of the spectra of free radicals 
with a rotating methyl or -C^ group may be predicted by a computer 
program based on a mathematical model developed by Anderson [ll"] and 
modified by Sack [12"]. 

The phenomenon of narrowing of magnetic resonance line breadths 
takes two forms; motional narrowing and exchange narrowing. The mecha- 
nism is in principle the same in both motional and exchange narrowing.. 

The model herein developed assumes that the precessing moments 
give rise to a radiated electromagnetic wave which is undergoing 
frequency modulation, because the magnetic interactions act to change 
the frequency of the precession. The frequency modulation is chang- 
ing in a random way in time due to the effect of the non-magnetic 
motions on the magnetic interactions. 

It can be shown in general that the spectrum of the radiation from 
any quantum-mechanical system is given by 

* ioot 2 

(1) I(oo) = Trace |\ M-(t)e dt | where h(t) is the radiating 

-00 

dipole moment matrix in the equation: 

(2) ih ^ = H|jl - p-H where H is the Hamiltonian. 

In any problem involving the type of narrowing dealt with here the 
Hamiltonian may be split up into three parts: 

H = H + H + H 
o p m 
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These three parts are the unperturbed Hamiltonian, the perturbing 
Hamiltonian- -generally just the dipolar interactions between the moments- 
which does not commute with H and thus can change the frequencies radi- 
ated by the system over some known range, and the.motional Hamiltonian 

whose characteristic is that it commutes with both H and u, and thus 

o 

can have no direct effect upon the radiation emitted or absorbed by the 
system, but does not commute with H and thus, by the relation 



(3) 



iftH = [H,H ] = [H ,H ] + [H ,H ] 
p p o p J m p J 



brackets denote the commutator can cause a time dependence of H . It 

P 

is this time dependence which narrows out the line broadening which H^ 
otherwise would cause. 

The random frequency modulation picture of the narrowing process 
can be derived from the following physical assumptions: 

(a) equation (1) 

(b) equation (3) 

(c) H - p,H =0 

m m 

(d) H H - H H =0 

mo o m 

(e) and, in addition, the assumption that H is small enough that 

P 

it has not important matrix elements connecting different states 
and E of the unperturbed energy H q . 

In any real substance, assumption (e) will not be rigorously true, but 

may be looked upon as an approximation which, if carefully handled, 

leads to good quantitative results. Under the above assumptions, the 

spectral intensity (1) will be computed. Consider one spectral line so 

that only the element of p.(t) which connects two unperturbed levels of 

H E. and E. need be considered. Then 
o 1 j 
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( 4 ) 



I(w) = 

i-j 



00 



-00 



, N -itut , i 2 
(j, . . (t)e dt 
ij 



Since that part of the time dependence of p» (t) ^ ue to unperturbed 

Hamiltonian is known, the following transformation may be performed: 



i , . , . ( 0 ) . 

M- . . = p. . .(t) exp(iut). . t) 

ij ij v v ij 



where 



Now 



, (0) = E. - E . 

hu). . ' l j 

ij 



+ H , p] = H p/ - p.' H 
dt m P P P 



because of assumption (a) and the fact that the transformation has re- 
moved the terms of the time equation (2) for p which depend on H q . By 

f 

use of assumption (e) , |Jk (t) ma Y t> e calculated as follows: 



ifili . ! 

ij 



= tv' - Vu 



(H ) . .Ub! . - p,! . (H ) . . 
p ll^lj lj p jj 



where 



Then 



\x ! . = - iAto. . ( t) |jl ! . (t) 
ij iJ iJ 



H (t). . - H (t) . . 

Aco . . ( t ) = - P - ■ — P- — II 

IJ h 



p,!.(t) = p,. exp(-i\ Au). . (t')dt ') 



ij 



ij 



ij 



This is inserted into the intensity formula (4) which gives, for the 
lineshape 

(5) I_(co) ~ exp(-i (cu - co.J 0 ^)t ” Aco^. ^ ( t ' )dt ' ) | 
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Aco (t) is a random function of time; its value at any time depends on 

ij 

the values of the diagonal elements of H at that time, and these will 

P 

change in a random way at a rate controlled by the motional Hamiltonian, 

as is seen from the time equation (3) (there is no time dependence due 

to H because these are diagonal elements) : 
o 



iftH = [H t H 1 

P m p 

The reason for assuming randomness is that in the important cases 

the effect of H back on H can be neglected, so that the motions 
p m 

embodied in H appear to the magnetic quantities to be uncorrelated, 
m 

Equation (5) and the idea of Aco„(t) as a random function of time are 

the basic ideas of the random frequency modulation model of narrowing 

phenomena. The term Ato_(t) is simply a construct, a single function 

which is meant to embody all of the broadening effects, diagonal as 

well as off diagonal. In order to find the line shipe, the assumption 

is made that Auo..(t) can be treated as a Markoffian random function, 
ij 

Equation (5) may be simplified in the following way: Redefine 

’ making U) the frequency deviation from the line 
center. Introduce the phase deviation T] ( t ) , which is the integral of 
the frequency deviation: 



cu as to = co - co. . 

ij 



T|(t) = \ Aco. . (t 1 ) dt ' 



0 



ij 



Then (5) becomes 



I (to) 



-S 



-TUt)) dtl 2 

-00 



Because of convenience, the correlation function form of the Fourier 
integral will be used. It can be shown that 
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or, 

( 6 ) 



IW = f e^dxf - 1 K t ’»dt 

-CO -CO 

I (to) = ^ e'*’ U)T cp (t ) cIt 



9(t) is called the correlation (or autocorrelation) function and can 
be written as 



( 7 ) 



, . , i(TKt-h-)-TKt). 

cp(t) - (e ) ave. over t 






t+T 



= <exp(i\ Ato (t ')dt ')> ave over t 



J t iJ 



This represents the averaged "memory" of the function at the time 

t+T for what its state was at a time T earlier. The problem is to 
find cp ( t ) under the assumption that A^k .(t) is Markoffian. 

A Markoffian function is a function for which the probability of 
a given value f 2 a t the time t, if the value was f at t - At, is inde- 
pendent of the value of the function at any earlier time than t - A t. 
Thus the probability depends only on the value of the function at the 
earlier time and not on its past history. The notation used is 

Probability = VKfJf^At) 

If the function is not to vary in a senselessly rapid fashion it is 
known that W must be proportional to At for small At unless f ^ 
Therefore, W may be written as 



(8) W(/ 1 |f 2 ,At) = 6 (f 1 ,f 2 ) + Ti(f I ,/ 2 ) At 

where rr is the probability per unit time of a transition from f to f 
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To solve the stated problem, a matrix method of solution is used 
in which the probabilities W ( f , /^ ) are considered as two dimensional 
matrices connecting the sets of indices to f ^ • The assumption is 
made that Aid has a discrete set of values. 

T 

The average in (7) may be written in another way. The interval T 

T 

is divided into n equally spaced steps, such that — is very small com- 
pared with the rate of change of Acd. Then 

n - t+T 

(9) exp iE — Au)(t+m^) = exp i\ Aoa ( t ' ) d t ' 

m=l J*- 



The probability that Acd takes on a certain set of values at the various 
instants in (9) is 



(10) W (Ail) ) W(Au> 



Alw 2 ’n^ xW ^ Atu 2 



Aco 3 ,J)x 



xW (Acd 



n- 1 



Acd , r ) 
n n 



cp (t ) may be found by summing the product of (9) and (10) over all the 1 
possible combinations of Aid's since these are respectively the value 
and the probability of a given averagand 



n 



n 



cp(r) = E E E *•* (exp(i tt E Aid )) xW (Acd ) tt W(Acd 

»»» 11 i m LL„ m- L 



Aid 



Acd Au^Acd^ 



m=l 



m=2 



m 



(ID 



n-1 

= E E • * * E tt (exp(i~ Acd ) x W(Acd I Acd ,— )) exp (iAoj — ) 

n m m ' m4 - 1 n r ' r> 

Aid^Acd 2 m=l 



m 1 m+l’n' 



nn 



The central product has the form of a matrix product. It will have this 
form exactly if the diagonal matrix to is defined as follows: 

Acd _ = 6 (Aid Aid . ) x Ato . 

=ij r j i 

The probability matrix is 

W (£) = W(Acd.|Acd., J) 
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3 H 



W^(Aw) is a row vector in Aw space, and J. (Aw), the vector all of whose 

components are unity is a column vector. 

T T 

Expanding exp(iAw^ — ) and neglecting Aw^ — with respect to unity, 

(11) is written as a matrix product multiplied on left and right by 
row and column vectors respectively: 



(12) 9(t) = W 1 • (exp(i — Aw) W(— ) ) n ^ * 1 

J- n = — n — 



Since ~ Aw can be made as small as is necessary, it is permissible to 
let 



T T 

exp(i — Aw) = 1 + i — Auo 

n — n = 



while from (7) the matrix W is rewritten in terms of matrix tt 

w(h) - i + 5 S 

This gives for (12) 



(exp(i J Aw) W(^)) n 1 = (1 + ^ (i|w + rr) ) n 1 



which if n goes to infinity is 
exp(T(i Aw + tt) ) 
so that (11) becomes 



cp ( t ) = • exp(T (iAw + tt) ) • _1 



or 



cp(T) 



- Z/ W.(Aco.) iexp(T(iAw + tt) V Aw Aw 
i,j _i 1 1 = ='J i j 



(13) = W exp[T(iQ + tt ) 1 ’ J; 

where the n components of the vector W are proportional to the 



87 



occupation probabilities of the sites in equilibrium, _I is a vector 
with all components equal to unity, Q is a diagonal matrix with ele- 
ments 0). and n is a matrix with elements 
J 



tt .. = p 
jk jk 



rr . . = -Sp,, 
JJ k J k 



j # k p = probability 



For cp(T) given in (13), equation (6) may be integrated as follows: 

\ exp( -ici)T)cp(T)dT = 

J o 

T 

W * [exp|r(TT + iQ - iioE^)^-]^ (tt + iQ - iiuE^) ^ * _L 

where is the unit matrix of order n. For T-*« the integral converges 
only if the real parts of all the eigenvalues of i Q + tt are non- 
positive. From physical arguments it is clear that the integral must 
converge, so that 



(14) I(to) ~ 2 Real W • (iuoE n - iQ - tt) • 1 



The calculation as to how the shape of a composite spectral line 
depends on the transition probabilities between the sites is thus 
reduced to the inversion of a matrix with complex elements. A computer 
program to perform the calculation in (14) has been written and begins 
on the next page. 
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"inn mo nrnr>o~>rn"> nnonnnnnn n noo n'o'o'noooocnmonooooooooo^o 



THIS PROGRAM CALCULATES THE I N T F MS I T Y 0 IS T P I RUT ION OF A 
rpMFpcn? SPECTRAL L I N c WITH EXCHANGE NAPROWING. 

R F F F p r nc F S FOR THF T HEOR Y APF 

J PHYS c 0C JAPAN, 9, 98° 

M n LFf I It A o PHYSICS ,1 ,1*3 



WRIT T PN PY RORFRT G. G T c C H APRIL 1969 



l^T C A P D RF T I POR PSP 2 POP NMR/NO. OF PLOT c ppo 
GPAPH/NO. OF GPAPHS/NO. OF POINTS TO PE PLOTTED/ 
width 0 c FREQUENCY SCALP IN INCHF~/DT WENS ION OF 
MATPIX/NO. OF NON7FR0 f) P p D I A GONA L EL F M EN T S ABOVE 
ni ACONf L/LINF WIDTH IN INCHESITHE RECIPROCAL OF 
T2)/THF HEIGHT, IN INCHES, OF thf HIGHEST ABSORPTION 
r FAK 11,12,12, 13, 312, 2F1 0.8 

2ND CARD VFCTOF OF INTENSITIES 6F10.5 

3 ° D CARP D T AGONAL MATRIX nr FREQUENCIES ce RESONANCE 
LI MP S . MATRIX SLPMENT AMO POSITION OF LINE IN 
I NC HF c 6(211, FI 0.5) 

4 T H CAPO THE TRANSITION PROBABILITIES TRANSFORMED Tn 
"HE T NC HF S c C A L F AC IS THE other data. the FIRST 
VAIUF CP Thf tr c NS I t T ON PROBABILITY FF D IN WILL 
HFTCPMINP THC HEIGHTS OP TH c AB c OR°T [pn FOP 

THE PF ST CF tHE PLCJR 6F10.5 

t TH CAPO N C N 7 P P P MA T P I X ELEMENTS ABOV c DIAGONAL 
7CI 1 

6 T H CADO LABELS FOR EACH PLOT 15 L ARFL C ON EACH 
Ci on 15A4 

7TH CAPO ADDITIONAL INFORMATION. TY °F IN A c MANY 
Cts R D S A c THFPF «pr CR A D HS Afa 



IMP Lie TT CQMCLF X*lfe (U,0,V,X) 

IMPLICIT pcAL*°(i-H,n,c,P-T,M,Y-7) 

Q F A L #4 XPEQ,XINTFN,YTNTEN 
R PA L LBI/4H / 

PFAI*R I TL ( 1 2 ) / * BOX 24, TOLL r P 

l • / 

DIME M Slew x< 20,2 0) , Omega (20,20) ,UO M EGA( 2C , 20 ) , 
lDDfiPf 20, 20 ,UPROR (2C ,20 , 1 Ml ( 2 0 , 20 ) , U f 20 , 20 ) , 

1FRFQ( l 00 r ) ,W(20) ,UW(20) ,UFREon OOO ) ,VELEM(20 ) , 

17 INTO N( 10C0) ,HI mten( 1000) , I A R ( f 0 ) ,BRROB( c O) , 

1XP PQ ( l OCO) , X I N T E N ( 1000) ,YI NTENf 1000 ) 

SF T matrices tc 7EP0 

DO 9 1=1,20 

on c j = 1 ,20 

OMEGA ( I , J) =0.0 

°POP ( I , J) =0. 0 

UOMFOA ( I , J)=(0. 0,0.0) 

U°P OB ( I ,J)=(G. 0,0.0) 

U( I , J )=(C. 0,0.0) 

9 UU( I , J) = ( 0.0, 0.0) 

E°=l. OE-BO 

S T A ° T TP* PEAD 7N DATA 

PF*D lCO,KKL,LLL,NNN,NUMPTS t I H I GH , N , MMM , C 0N C T , R CALF 
READ 110,(W( I) ,1=1 ,N) 

F R E 0 ( 1 ) = 0 . 0 
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non ~> -> ") im 



MM] = M ( j w c T c - 1 
iH T G»-*= THIGH 
ANMPTc = miimot': 

7 I MT =/ u]r,|j/HMrTc 
no ] T - ] f MM] 

1 cproi H|)--orr)|i |,7TMr 

D c .' n I ;> , ( « r , j ,m r,/u',ni,^i,M 
I. LMf\'=l LI *N\'v 

o-AH 1 1C, C (BFP^P (T ) ) , I =1 »l.LMM) 
nn l? c T =1 ♦ L t. N N 

12 = RCpnR ( r j ; p^op ( t ) f^ % pq 31 a 6 
MM= pj tMK’M 

OF'P 10 c , < t T£P( I ) ) ,1=1 ,mm) 

P r fK» t-t tlF fico, 12 ) 

P'f | *0 P I A N* ( 1 2 ) / ’ 

l » / 

P r -t I ■ ' nf I ( IOC ) 

RFA P 12r ,( L a RF! (?) ,1=1 f L l. N M ) 
nn 2 r r T=i f M\N 

300 P^n 121t(? T ? T l p (T * J ) » J= 7 » 1 2 ) 

CMP P r / ^ ]v PATA, START 0/ L CMLA ‘ r T n MR 

PH 1 r ^ 1 JT =1 ,MVV 

on 1C n] t |r] f i | l 

czr up cprp/RTLr r Y MATRIX 

nn cc T = 1 » mm f 2 
11=1 + 1 
i * = I /. p ( n 
TP=TA R ( T 1 ) 

r S OPHBI I' » ! P ) =-R P ROP ( ( JT-1) *Ll L+T J) 
no cc t = i t M 
P° j = 1 » N 
cc p R o n ( j , ] ) =oo rp ( I t J) 

0° 22 T = 1 » N 
2 ? ^PP° ( T t J )=C.C 

no ^6 t - i ♦ n 

7=0.0 

no r 7 J=l t N 

” 7 =7 + r, R r R ( I , J) 

c F soflRf T ,? ) = -7+C0N c T/6. 2R31°s 

CCT !|C P ? f 0 C M A L Toampjttcim F'O'QIJRN'CY matMX AN'O 
pnMpTf.r fF-"i j-tf matc«IC c S t 0 RTS i JIT ANT m/.tpjx 

nn ac T = 1 , N 
on . 1 = 1 » n 
u^Mf ct ( r ♦ j)=cmfoa ( t ,. n 

( i , j)=uomc G a d f j) *(c.o,-i .0) 

IJPROP ( T ♦ J) =oonp (I f j) 

JIPPPP ( T t J)=MDPPP(I f J)*<-1 .c f 0.C ) 

' C UU( T ♦ J)=U0* c G* ( T , J ) + 1 1 c p n fj ( t f )) 

no LOOP UH T r H r./LCULATrS IVTFNSITT C ^ 

on 1CCC M = i f M|MPTS 
on 11 1 = 1 t N 
0° 11 J=1,N 
11 U( I tJ)=UIJ( T , J) 

IJFPCP( N)=FRR0( M) 

llFPCOf m)=|jprfQ( m) *<C. C»1 .0) 

cpcoilRMCY H=mpmt A POF 0 TO imaginary °ART OP p T AO.OM AL 
FI .c Mr NT c 

on 2 T=1,N 

2 U( I ♦ T ) =(K I t T ) + 1 ic cp q ( m ) 

MATRIX T^!V/P D FI^N 
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non non non ~>or> ono 



TALL GA('SF 4 (N!,FP,U,X,KFQ) 

VFC Tno PF INTPNSI TTF C *I NV c PTFO MATRIX 

DP 6 J = l f N 
0=( C.Ct 0.0) 
nn *> k = i,n 
UW( K)=M K) 

V M ULT=UW ( K) *X( K , J) 

5 o=o+vmiilt 

6 VFL c M(J)=0 

RE SUL T^UMI T V c CTrp 
n = ( C.C,C.O) 

nn 7 I = 1 , N 

7 q=o + VFI.fm( t ) 

Ti'k’f.r pfal PART P c THI C R c 5 U t.T 
p=o 

INTENSITY TS t MCF THF RFAl PART CF THIS RFS'JL' p 
1000 7 INTFMf d= 2 . *S 

C ND CALCULATIONS. STAPT COMPUTING PRAPHS 
I F ( KKL.FO. n OP TO 450 

NMD CCfTrpM 

on 5 C C I=1,M 
X° F Q f I )=PPCf)( I ) 

5 C 0 X! MT F N( I ) =7 1 NTFNf I ) 

I c f JI .GT. 1 ) GO T P 

I f f i j .OT. i ) nn tp, o< 5 P 
T np=C . 0 

or 444 . T = l t NljMPTC 

I F f TPF.OT.XINTFNII )) on TQ 444 
TPP = X I N T F N ( T ) 

4*4 CPNTIWUF 

qcq p n 454 I=1,NUMPT$ 

454 XIMTFNC I )=X! NTFN( n*SCALF/Tpo 
IFf TJ.FO. DOC T P 765 
OP 766 1 = 7,12 
7*6 T T L C I ) =P L A NK ( T ) 

OP Tp 767 
765 n n 7 *o 1 = 7,] 2 
^69 I T L( T ) =ITITL.F (JI ,1 ) 

7f7 CONTI NLF 

t.BL=t.AFFL( ( JI -1 ) *11 L+ I J ) 

TF( LLL.ro. l) on TO 204 
IFf TJ.FO. 1) on tp 201 
I c ( I J.FO. ILL) 00 TO 202 

CALL OPAW( NU"DT.S,XINTFN,XRFQ,? ,0,LBL,ITL, 1 .C, 1.0, 0,0, 

12. 2, °, THIGH, 0 , LA ST) 

00 TP 1001 

2d CALL DRAW(NUMPTS,XINTFN,XRFQ,l,O,LRL,ITL,1.0, 1.0, 0,0, 

12. 2, °, THIGH, 0, LAST) 

OP TP 1 C 01 

204 CALL nPAW(NI)^P T S,XINTEN,XRFO,o, 0 ,LBL,ITL,l . 0 , l.C, 0 ,C, 
12 , 2 , 9 , THIGH, C, LAST) 

on tc 1 cm 

202 CALL DPAW(NU m P t . c ,XTNTFM,XP c 0 , 3 , 0 ,LBL,I t L, 1 . 0 , 1 . 0 , 0 , 0 , 

12. 2, f, THIGH, 0, LA ST) 

On T P 1 CC 1 

F S° SFCTlfN 

450 DP 9 05 1=1, WMl 
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C CMC'JLATpc NUMERICAL DERIVATIVE 
C 

30^ H I N T r N ( T ) = ( 7INTFM( 1 + 1 )-ZIN tc M(T))/7INC 
nn *> c l T = 1 ,M 
yorof I )=rocQ(i ) 

S r l Y T A ,T r-N*( T )=HI MFMI I ) 

iPCJl.r.T.D nn th qqc 
! f ( I j ,r,T, i ) nn t c qqo 
t ° p = c . r 

HO 445 

i r ( tod .gt.arr( yintfn( 1 ) ) ) nn T p 445 

T^PrARSI YINTCM( I ) ) 

44P rnNIJT Ml IF 

ICQ nn 4CC T=i f ^Ml 

4 C 5 Y T N T T M ( T )=YTNTFM(I ) *5 C & LF / ( TOP*? . ) +4 . 5 
I F ( IJ.PO. 1 )G0 T G 86F 

nr 86* 1 = 7,12 

* DL( I )=PDNK(I ) 
no TO 0f7 
a A* no P/: q [=7,1?. 

8*8 TTU T )=TTITLF'( JI t I ) 

8A7 CON'TINLT 

LA|.=LAR r L( (.11-1 )*LLL+I J) 

I F ( LI L.FQ. 1 ) on Tr 304 
J r ( IJ.ro. 1 ) nn TO 30 1 
T f ( i j.fo.lld nn to 302 

CD L n P/ t W ( N Mi,YIMTFM,XRFQ,2,0,LBL,ITL t l .0,1. 0,0, 0,2, 2, 
P, IHTGH,0,Lt r *T) 
no to i C01 

301 CALL D D A M N M 1 , Y ! NTE N , X 0 F 0 , 1 , 0 , L R L » I T L » 1 . C , 1 . ") » 0 » 0 , 2 * 2 » 
l *, I HI OH, C, LA ST) 

no ' r n icci 

30* CALL f)C*W( N’Ml ,YINTFN,X°EO,0,0 ,LRL, ITL,1.C,1.0,D,0,2,2, 
1 F , I h I OH , C , LA ST) 

GO T P 1001 

302 CM L DPAW(NM1 ,YINTON,XFF0,3,3,LRL, ITL,1 .0, 1.0,0,0,2t2, 
1P,IHIGH,C,LA r T) 

loci cnMTTNijr 

100 FQPm*t( 11,12,12,13*312 »2F 10. *>) 

102 FORMA T ( 6 ( 2 1 1 ,F 1 0. 5 ) ) 

ICS forma T( 7C( ID) 

110 FP.PMAT(6F10. 5) 

120 Forma t( 1 c A 4 ) 

121 FOPMA T(6AP) 

PToo 

ENir 



SURRGUTINfF GAUSS4( N^o ,a ,x,rcp) 

IMPLICIT C0mplFX*D ( 4-n,F-H,0-Y) 
niMCNJSTCM 4 ( 20 , 20 ) ,X( 20 , 20 ) 
on 1 1 = 1 , a.' 
on 1 j = 1 , n 

1 X( ! ,J )=( 0 . 0 , 0 . 0 ) 
no ? K = 1 , N 

2 X( K ) = ( l . 0 , 0 . 0 ) 

1 C DO 34 L = 1 , N 

K P = C 
7 = 0 .c 

on 12 k=l,n 

I F ( Z-CDABMA (K,L) ) ) 11 , 12 ,12 

11 7 =COAR c (A(K,L)) 

« D =K 

12 CONTI MUO 
IF(L-KP) 13 , 20,20 

13 DO 14 J = L , N 
S=A(L, J) 

A ( L , J )=/• ( Ko, J) 

14 A(KP,J)=$ 
no 1 * j= 1 , n 
S«X(L, J) 

X( L , J ) = X( K p t J ) 
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1*5 X(KP,J) = c 

20 TF(rOARC(ML»L) )-F D )FQ, 50,30 

30 TP< L-N)3l ,34,34 

31 LP1=L+1 

nn 3 A k = t o i f Ni 

TF(CP£ A ( K ,L) 1 >32,36,32 

32 RA Tin = ^ ( K , L ) /A ( l , L ) 

DO 33 J = L c> 1 , N 

33 A ( K , J )=* ( K, J )-p.ft-rin*A ( L j) 
no 3* J = 1 , N 

3 F X( K, J) = X( K, l)-P/TIC*X(l , J) 

36 rONTINUF 

34 fOMTl \u;r 

4C no 43 I =1 ,N 
TI=M+1-I 
no 43 J = 1 , N 
"' = { 0.0, C.0) 

IF ( I I-NK1 ,43,' 3 
41 ITPl=!I+l 

on 4? K=TI P 1,M 

4 2 c = ( 1 1 ,K) *X( K , J) 

43 Xf ! T , J)=( X( II , J)-5) /MI I ,1 !) 
«fo = \ 

RFTIjon 

K p P=2 

R F T | |p M 

cfin 
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